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Further one-loop results in 0(a) improved lattice QCD0 
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Using the Schrodinger functional we have computed a variety of renormalized on-shell correlation functions to 
one-loop order of perturbation theory. By studying their approach to the continuum limit we have determined 
the O(a) counterterms needed to improve the quark mass and a number of isovector quark bilinear operators. 



1. INTRODUCTION 

Recent work by the ALPHA collaboration has 
focused on non-perturbative renormalization and 
on-shell 0(a) improvement of lattice QCD with 
Wilson quarks [1-7]. Various improvement coef- 
ficients could be determined as functions of the 
bare coupling g 0l including the coefficient of the 
Sheikholeslami-Wohlert term in the lattice ac- 
tion g. 

While a non-perturbative determination of the 
improvement coefficients is clearly preferable per- 
turbative estimates are nevertheless useful. First 
of all, any non-perturbative determination should 
establish contact with perturbation theory at suf- 
ficiently small values of the bare coupling con- 
stant. This provides a non-trivial check for the 
chosen strategy and criteria to assess the quality 
of a given improvement condition. Secondly, for 
those coefficients which account for lattice effects 
due to non-zero quark masses, perturbative es- 
timates may indeed be satisfactory provided the 
quark masses are small when measured in lattice 
units. 

In this contribution we present our one-loop 
results for the on-shell 0(a) improved isovector 
composite operators which are bilinear in the 
quark fields. The computational strategy and 
most of the results have already been published 
in refs. PJlCj]- In addition we here also include 
the results for the improved isovector tensor and 
scalar densities. 



2. DEFINITIONS 

We consider lattice QCD with Nf > 2 degen- 
erate quark flavours of bare mass mo and shall 
assume that the action has already been on-shell 
improved We are interested in the improve- 
ment of the following isovector operators, 



P a (x) 
S a (x) 



^{x)\r a ^{x), 
i4'(x)a^^T a ifj(x). 



(1) 
(2) 
(3) 
(4) 
(5) 



Here r a (a — 1, 2, 3) arc the usual Pauli matrices 
acting in any two-flavour subspace and our con- 
ventions for the Dirac matrices are as in ref. ||. 

In a mass independent renormalization scheme 
the renormalized (at renormalization scale /i) and 
0(a) improved counterparts of the above fields all 
take the form (A = V, A, P, S, T) [§ 
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A R = Z x (<?o, an) [l + 6 x (3o) am c,. 
where Ai stands for 
{Vi)l = V; + cv(g 2 )ad u T^, 
(Aj)l = Al + c K {gl)ad^P a , 
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(2i)J„ 



•or(sg)a(W -d v V?), 



(6) 

(7) 
(8) 
(9) 



and otherwise X\ — A. Here denotes the sym- 
metric lattice derivative and the parameter go is 
connected to the bare coupling go through 
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90 = 00 [1 + b g(9o) am q] 



(10) 
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where m q = mo — m c and m c is the critical bare 
quark mass. Similarly one defines 



m q = niq [1 + 6 m (go)aw q ] , 



(11) 



and the renormalized 0(a) improved coupling 
and quark mass are then related to these param- 
eters by 0, 



9k = <7o^gG?oW)> 
tor = rh q Z m (gQ,a^). 



(12) 
(13) 



The improvement coefficients & g , o m , 6x can be ex- 
panded in perturbation theory, 



& = 6(°)+&« 9o 2 + 0( 5o 4 ), 



(14) 



and an analogous expansion exists for the coeffi- 
cients ca, cv and ct- 

3. COMPUTATIONAL STRATEGY 

To determine the improvement coefficients we 
chose to compute a number of on-shell correla- 
tion functions derived from the Schrodinger func- 
tional (SF). The SF is the Euclidean functional 
integral for QCD on a finite space-time manifold 
where the (spatially periodic) quantum fields sat- 
isfy Dirichlct boundary conditions in the time di- 
rection (nj[l§]. For proper choice of the boundary 
conditions it can be shown that the lattice action 
has a unique absolute minimum p| . The sad- 
dle point expansion about this minimum is then 
straightforward (albeit technically involved), and 
zero modes do not appear. 

The gauge field boundary conditions imply that 
only global gauge transformations are allowed at 
the boundaries. Therefore, gauge invariant corre- 
lation functions can be defined where the quark 
and antiquark fields at the boundaries are sepa- 
rately projected onto their zero spatial momen- 
tum components. This is convenient because the 
perturbative expansion of such a correlation func- 
tion starts with tree diagrams. Furthermore, ex- 
actly the same correlation functions can be used 
in numerical simulations. 

In order to take the continuum limit in a fi- 
nite space-time volume one fixes the time extent 
T, the renormalized 0(a) improved quark mass 
and all other dimensionful parameter in units of 



L, the spatial extent of the space-time manifold. 
As a result 0(a) lattice artefacts always appear 
as a/L effects and can be identified by varying 
the lattice size. In each renormalized correlation 
function these effects are cancelled by an a priori 
different linear combination of 0(a) improvement 
coefficients. The finite volume provides a great 
flexibility here because many different renormal- 
ized correlation functions can be obtained by sim- 
ple changes of the boundary conditions. 

For completeness we mention that Dirichlct 
boundary conditions cause additional divergences 
and 0(a) artefacts localised at the boundaries. 
These can be absorbed by renormalizing the 
boundary quark and antiquark fields in the same 
way as the composite fields in eq. (^), and by in- 
cluding additional 0(a) boundary counterterms 
in the lattice action |ll] , |l2| ,|3[ . 

For details of our computational strategy and 
the definitions of most of the correlation functions 
the reader should consult refs. We have 

treated the case of the isovector tensor density in 
complete analogy to the improved vector current, 
and the computation of 6s involved a correlation 
function similar to /y of ref. ||, where the vec- 
tor current was replaced by the isovector scalar 
density. 
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Table 1. Improvement coefficients b 



4. RESULTS 

To one-loop order of perturbation theory we 
have carried out many consistency checks and 
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thus confirm the general framework of 0(a) im- 
provement as described in ref. j|. 

Numerically we obtain (C F = (N 2 - l)/2N for 
N colours), 

c A = -0.005680(2) x C F gl + 0(g%), (15) 
c v = -0.01225(1) xC FSo 2 + 0(^), (16) 
c T = 0.00896(1) xC FSo 2 + 0( 9o 4 ). (17) 

The coefficient ca has first been obtained in 
ref. M . and cy was given in ref. jl0| . The results 
and references for the ^-coefficients are collected 
in table 1. Note that to order gfi we find, within 
errors, 

b s = -2b m . (18) 

In fact, it can be shown that eq. ( |l8| ) is an exact 
identity in quenched QCD and furthermore the 
isoscalar scalar operator has the same ^-coefficient 
as the isovector Q. 

Comparison with the non-perturbative results 
for ca and c v shows that in these cases 
perturbative estimates are not accurate at large 
values of the bare coupling constant. In the case 
of by we can compare with the non-perturbative 
result of ref. || . In figure 1 we see that the non- 
perturbative values (represented by the solid line) 
are quite a bit higher than the perturbative es- 
timate (dotted line), even when Parisi's boosted 
coupling Q is used (crosses). However, using the 
boosted perturbative one-loop estimate at go = 1 
an error of only 1 per cent is induced in the nor- 
malisation factor 1 + byam q [cf . eq. @] , provided 
am q < 0.05. 

It is of course not clear whether similar con- 
clusions can be drawn for the other ^-coefficients. 
Further non-perturbative results would obviously 
be welcome and some progress in this direction 
has been reported at this conference ]lq , p7t . 
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